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Roles of multipole degrees of freedom in multiband superconductors are investigated in a 
case of impurities whose low-lying states consist of singlet ground and triplet excited states, 
which is related to the experimental fact that the transition temperature T c is increased by Pr 
substitution for La in LaOs4Sbi2. The most important contribution to the T c increase comes 
from the inelastic interband scattering of electrons coupled to quadrupole or octupole moments 
of impurities. It is found that a magnetic field modifies an effective pairing interaction and the 
scattering anisotropy appears in the field-orientation dependence of the upper critical field 
H C 2 in the vicinity of T c , although a uniaxial anisotropic field is required for experimental 
detection. This would be proof that the Pr internal degrees of freedom are relevant to the 
stability of superconductivity in (Lai_ a; Pr :! ;)Os4Sbi2. 
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. 1. Introduction 



It is common knowledge that magnetic impurities in- 
duce pair breaking and a decrease in superconducting 
transition temperature T c . In a multiband case, however, 
there is room to reconsider such a conventional under- 
standing. 1, 2 ' Electron scattering by impurities occurs not 
only within one band (intraband) but also between dif- 
ferent bands (interband). In addition, variations in order 
parameters can be considered for superconducting bands. 
In a two-band case, for instance, the s±-wave state is 
characterized by the sign- reversing order parameters. 3 ' 4 ' 
According to recent studies on the s±-wave state, in- 
terband scattering does not contribute to T c suppression 
for magnetic impurities. 5 ' 6 ' This behavior resembles that 
of intraband scattering by nonmagnetic impurities in a 
single-band s-wave superconductor. Furthermore, T c can 
be increased by inelastic interband scattering if the im- 
purities have such internal degrees of freedom as crystal- 
field split energy levels. 7 ' The idea was originally pro- 
posed by Fulde et al. for nonmagnetic impurities in a 
single-band superconductor; 8 ' however, it has been con- 
sidered that magnetic impurities always induce T c sup- 
pression. We have to check carefully the potential roles 
of impurities with orbital degrees of freedom, e.g., mul- 
tipoles in /-electron systems. 

Since the discovery of the heavy fermion supercon- 
ductor PrOs4Sbi2 (T c = 1.85 K), 9 ' much attention has 
been attracted by the following unique superconducting 
properties: 10 ' (1) Change in gap symmetry in a mag- 
netic field. This was first observed as the field orientation 
dependence of oscillating patterns on thermal conductiv- 
ity, 11 ' although no symptom of cubic symmetry breaking 
was found in a recent angle-resolved specific heat exper- 
iment. 12 ' (2) Broken time reversal symmetry below T c . 
An internal magnetic field emerges spontaneously but is 
extremely small, which was detected by muon spin relax- 
ation measurement. 13 ' The identification of the pairing 



parity is under debate. (3) Robustness against substitu- 
tion effects. Substituting La for Pr leads to the grad- 
ual decrease in T c . 14 ' Surprisingly, the x dependence of 
T c in (Pri_ ;E La 2 ;)Os4Sbi2 is smoothly connected to the 
T c (= 0.74 K) of LaOs4Sbi2 that exhibits the conven- 
tional s-wave property. 14 ' 15 ' In another experiment on 
Pr(Osi_ :E Ru a ;)4Sbi2, T c gradually decreases until x ~ 0.6 
and then shows a monotonic increase up to the T c (= 1.3 
K) of PrRu4Sbi2, which is also an s-wave superconduc- 
tor. 16 ' It has been considered that these unconventional 
properties are associated with the Pr 4/-electron behav- 
ior, although this association is still mysterious. 

It is also important to note that PrOs4Sb 12 exhibits 
the quadrupole ordering in a higher magnetic field where 
the superconductivity disappears. 10 ' 17 ' 18 ' This ordering 
phase diagram is successfully explained by the localized 
Pr 4/ 2 quasi-quartet model. 19,20 ' In fact, it has been 
established by the recent experiments that the Pr low- 

(2) 

lying states consist of T± singlet ground and T\ ' excited 
triplet states well separated from the higher crystal-field 
energy levels. 21,22 ' The r^ 2 -* representation is for the 
point group with no fourfold symmetry axis. 23 ' This wave 
function is expressed by a combination of T4 and T§ wave 
functions in the Oh point group, where the latter is 
more dominant. Owing to the small singlet-triplet en- 
ergy splitting ~ 8 K, it is expected that the quadrupole 
degrees of freedom will play a key role in the unconven- 
tional superconductivity as well as in the field-induced 
ordering. However, there is no reason to deny any con- 
tribution from other multipole degrees of freedom in the 
Pr singlet-triplet configuration. 

To find a clue to understand the multipole contribu- 
tion to the superconductivity, we focus on the Pr impu- 
rity effect on the LaOs4Sbi2 superconductor. As men- 
tioned above, T c is increased by Pr substitution for La in 
LaOs4Sbi2- In addition, multiband properties are indi- 
cated by thermal-transport measurement under a mag- 
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netic field in PrOs4Sbi2 24 ^ and by nuclear quadrupolc 
resonance measurement in (Lai-^Pr^C^Sb^. 15 ) Sup- 
pose that LaOs4Sbi2 is a single-band superconductor 
and that the inelastic scattering by the Pr impurities 
contributes to the T c increase, the most probable ori- 
gin of the T c increase is a quadrupolar scattering effect 
that leads to an effective pairing interaction that stabi- 
lizes the superconductivity, by analogy with the optical 
phonon-mediated pairing. 25 ) If a multiband picture is ap- 
plicable, there are two possibilities: For the s ++ -wave 
state with the same sign order parameters, interband 
nonmagnetic (quadrupolar) scattering is also favorable 
for the T c increase by the inelastic impurity scattering, 
while interband magnetic (octupolar) scattering can in- 
crease T c for the s±-wave state. The relevance of the 
multiband scenario is closely connected to the character- 
istic structure of skutterudites. In PrOs4Sbi2, each Pr 
ion is located at the center of the Sbi2 icosahedron cage. 
The most important point is the local hybridization of 
the Pr 4/-electron states with conduction bands via the 
Sb 12 molecular orbitals denoted by the a u and t u point- 
group symmetries. 26 ' 27 - ) If the a u -t u orbital exchange is 
the most relevant for the Pr impurity scattering, the mul- 
tipolar coupling in the Pr singlet-triplet configuration 
contributes to the T c increase. In our previous work, we 
proposed a possibility of T c increase by magnetic impu- 
rity scattering in the s ++ -wave state erroneously instead 
of the s±-wave state. 7, 28 - ) We will give a correct descrip- 
tion of the impurity scattering effect on the two-band 
superconductivity in this case. 

In this paper, we discuss the multiorbital scattering 
effect on the multiband superconductivity that reflects 
the local orbital symmetry. In the Pr singlet-triplet con- 
figuration, electrons are coupled to the quadrupoles ex- 
pressed as yz, zx, and xy or the octupoles expressed as 
x(y 2 — z 2 ), y(z 2 — x 2 ), and z(x 2 ~y 2 ) by analogy with the 
dipoles expressed as x, y, and z, respectively. 29 -* We focus 
on the magnetic octupolar scattering effect throughout 
the paper. There is a marked distinction between such 
multipoles and the spin when a magnetic field is intro- 
duced. For the orbital exchange scattering by a multipole 
moment, we find that its polarization changes as the field 
direction is rotated, while the spin exchange scattering is 
isotropic. If such orbital exchange scattering is an origin 
of T c increase, an effective pairing interaction is mod- 
ified by the scattering anisotropy under the field. It is 
expected that the anisotropy effect will be observed as 
the field orientation dependence of an upper critical field 
H C 2(T) line. This would be conclusive evidence indicat- 
ing that the orbital degrees of freedom definitely play a 
crucial role in the T c increase, which is closely connected 
to the unique Pr atomic structure in (Lai-^Pr^C^Sb]^. 

The paper is organized as follows. In §2, we show a typ- 
ical form of the inelastic magnetic scattering that cou- 
ples the Oh Ti singlet and triplet states. In §3, we 
explain a gap equation for T c increased or reduced by 
interband magnetic scattering as an impurity effect on 
the two-band s-wave superconductivity. The same for- 
mulation is applied straightforwardly to an interband 
nonmagnetic scattering case. The gap equation is mod- 
ified by including a magnetic field effect on impurities. 



In §4, this argument is extended to a case of single a u 
and threefold degenerate t u bands that reflects the cu- 
bic symmetry. Owing to the Zeeman splitting of the ex- 
cited triplet, T c depends on the field direction. Since the 
T c deviation is very small when the cubic symmetry is 
conserved, we consider a uniaxial anisotropy effect in the 
gap equation to elucidate the close correlation between 
the field-orientation-dependent T c and the anisotropy of 
multipolar scattering. Finally, conclusions are given in 
§5- 

2. Inelastic Scattering by Impurities 

First, we show a typical case of the inelastic electron 
scattering by magnetic impurities, keeping in mind the 
Pr 3+ states in PrOs4Sbi 2 . We consider that the low- lying 
states consist of the Ti singlet ground and T 5 triplet 
excited states in an Oh crystal field. They are expressed 
by 29 ) 



(1) 



|r 1 ) = ^(|4) + |-4)) + ^l|0), 



ol-l). 

2)), 




(2) 



3) 



1). 



where \M) (M — 4, 3, • • • , —4) is an eigenstate of J z for 
the J = 4 total angular momentum. The interchange 
of the singlet and triplet states is caused by the local 
orbital exchange of electrons via the hybridization of /- 
orbitals with conduction bands. Since each Pr ion is lo- 
cated at the center of the Sbi2 cage, the Sbi2 molecu- 
lar orbitals mediate the hybridization, which is the most 
pronounced feature of skutterudites. The a u (T2) molec- 
ular orbital mostly contributes to the main conduction 
band (a u band), hybridizing with the /-orbitals most 
strongly. 26 ^ The spin and orbitally coupled states of /- 
electrons are categorized as the T7 representation. We 
also consider here that the t u (T4) orbitals participate 
in the secondary conduction band (t u band), having a 
weaker hybridization with the Tg /-electron states (see 
Appendix A). 

Both quadrupole and octupole moments with the T§ 
symmetry are involved in the interchange of the Pr Ti 
and T 5 states. We focus on the inelastic octupolar scat- 
tering since it is unconventional that T c can be increased 
by magnetic correlations. In the present case, the oc- 
tupolar exchange scattering is accompanied by the a u -t u 
orbital exchange of local electrons, which is shown in 
Fig. 1. This effective exchange interaction is expressed 
by the following Hamiltonian W. cx = V.i + V.': 



(3) 



W = J s ^2 X] / dra \n a ~fn>6(r - Ry) 
R-y nn' 

x^(r)(l s ) 4>(r). (4) 

V / nn' 
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Fig. 1. Sketch of the inelastic multipolar exchange scattering 
with the a u -t u orbital exchange in the Pr singlet-triplet con- 
figuration. Here, the local t u and a u electrons participate in the 
different bands denoted by (+) and (— ), respectively. 



Table I. Effective pairing interaction type mediated by the inter- 
band magnetic or nonmagnetic (octupolar or quadrupolar, re- 
spectively, in the Oh Ti-Ps configuration) scattering impurities: 
A (attractive) or R (repulsive). It depends on the pairing type 
used. 



Scattering type 


Pairing type 


Interaction type 


magnetic 


s-j--)_-wave 


R 




s±-w&ve 


A 


nonmagnetic 


s.)--(--wave 


A 




s-l--wave 


R 



trons: 



H c = 



E 



drV>L(r)e(-iV)V>M<r( r ) 



dr 



The first term Hi is for the impurity states, where i? 7 
represents the position of the 7th impurity and at fn (a 7 „) 
is the pseudo-fermion creation (annihilation) operator for 
the nth impurity energy level 5 n at the 7th impurity 
site (n = 1,2,3, and 4 denote Ti, T 5 +, r 5 0, and r 5 — , 
respectively). 30 - ) In the second term W , Is represents 
the rVtype octupolar exchange scattering of conduction 
electrons with the coupling constant Jg (S denotes the 
magnetic scattering here) that is accompanied by the 
interchange among the nth and n'th energy levels: 



(is) 

\ f 1 



(Tz) nn , 4 + \ (T+) nn , L + \ (T_) w t +> 



(5) 



where T v and t n (n = z, ±) are octupole operators for the 
impurity states and conduction electrons, respectively. 
Their complete expressions are given in §3 of ref. 31. 
We use the two-band expression xjjH^ip for the octupolar 
exchange scattering with the a u -t u orbital exchange as 
described in Appendix A, where we have corrected the 
previous results. 7 ) 

3. Gap Equation for Two-Band Superconductiv- 
ity 

Before discussing the multipolar scattering effect on 
T c , we briefly review our previous study and give a cor- 
rect description for the the Pr-like impurities with the 
singlet-triplet configuration, focusing on how to derive a 
gap equation for the two-band s++-wave and s±-wave 
superconducting states. 7 ) A magnetic field effect, which 
modifies an effective pairing interaction, is also taken into 
account in the gap equation. In Table I, we summarize 
the effective interaction type, which is either attractive or 
repulsive, mediated by the interband magnetic (octupo- 
lar) scattering or nonmagnetic (quadrupolar) scattering 
impurities in each s-wave state. 

3. 1 Formulation 

To describe the superconductivity, we start from the 
following Hamiltonian for the two-band (fj, = ±) elec- 



(6) 

Here, we assume that both bands are identical and the or- 
der parameters have the same amplitude |A+| = |A_| = 
A for simplicity. The operator e(— iV) = — V 2 /(2m e ) — 
E-p expresses the kinetic energy measured from the Fermi 
energy Ep, where m e is the electron mass and h = 1. It 
is convenient to introduce the 8x8 matrix form of the 
thermal Green's function 



G(T,r,r') = -(T*(r,r)*T(r',0)) 



(7) 



with the eight-dimensional vectors \&(r) and \J^(r) for 
the two-band electrons defined as 



and 



/ Mr) \ 

Vvi( r ) 



t _ 



*1 



*L ), 



(8) 



m 1 - 



(r) = ( V^M V^(r) ^(r) ^(r) ) .(9) 
Their Heisenberg representations are written as 

Ut 
(10) 



* M (r, r) = e WT * AI (r)e"^, *t(r, r) = e HT ^Jr)e 



In the absence of impurity scattering, the unperturbed 
Green's function is Fourier-transformed to 



G (iwi, k) 



A 5 ' 



(11) 



where a a is the Pauli matrix for the spin space and 
p a is that for the particle-hole space (a = 1,2, and 
3 correspond to x, y, and z, respectively), and e& = 
k 2 /(2m c ) — E-p. The two-band superconductivity is ex- 
pressed with another Pauli matrix f Q in the band space 
as 



4±(l+T3) + ^(l-f 3 ) 



A I -!- r.) - — (I - r.,). (12) 

We consider the magnetic (octupolar) scattering effect 
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on T c when the impurities are distributed randomly in 
the two-band s-wave superconductor, which is expressed 
on the right-hand side of the following linearized gap 
equation: 7 ' 8 ) 



8T C 

7T 



T log 



T 

± c 



fu / A 
/a fu 



A_ 



(13) 



where T c (T c0 ) is the transition temperature in the pres- 
ence (absence) of impurities, /a and f u represent self- 
energies corresponding to the order parameter and Mat- 
subara frequency components, respectively. On the basis 
of the unperturbed Hamiltonian He + Wi in eqs. (6) and 
(3), the self-energy is obtained as 

1 1 



lWl 



1W 2 



x J; 



G (iw; + iwi - iw2, k) 



('«). 

V / n'n 



(14) 



in the second Born approximation by W in eq. (4) for the 
magnetic type of effective exchange scattering Is in 
eq. (5). Here, n{ mp is the impurity density, T is the tem- 
perature (fee = 1), and f2 represents the system volume. 
By analogy with the optical phonon case, the inelastic 
impurity scattering leads to an effective pairing interac- 
tion. After calculating the self-energy terms, 



Sa(Iwz) 



-Tr 



/S 2 <x 2 £s(iwi) , S w (iw ( ) = 



-TrSsK), 
(15) 



where 1 /8 is the normalization factor in the f <£> p ® a 
space (f for the band space, p for the particle-hole space 
and a for the spin space), we obtain 



8T c r 



-/a (a;) > 0, 



i S u (iwz) 7r 



8T c r 



T f u (x)>0, (16) 



where x represents the energy difference between the T\ 
singlet ground and T 5 triplet excited states as 

s r 5 - 5 Tl 



2T r 



(17) 



We define the lifetime due to the magnetic impurity scat- 
tering as 

^ = 27m imp JVo ^) ' ( 18 ) 

where No represents the density of electronic states at 
the Fermi energy. In the absence of a magnetic field, 
both /a{%) and f u {x) are independent of 9 and <p re- 
lated to the local hybridization defined in eq. (A- 10). In 
eq. (13), the matrix gives the higher eigenvalue f(x) = 
— Ia{x)+ fui{x) > for the higher T c in the s±-wave state 
(A + /A_ = —1). The explicit representation of f(x) is 
given in Appendix B. The s + +-wave state (A + /A_ = 1) 
is for the lower T c . In the same manner, for the nonmag- 



netic (quadrupolar) scattering in Appendix A, the gap 
equation is given by 



t w log — — 



c0 



Ah 

A_ 



fui "/A 
-/a fu 



A 4 
A_ 



(19) 

where is the corresponding lifetime, which leads to 
the T c increase in the s ++ -wave state (see Table I). 

3.2 Magnetic field effect 

Next, we consider a magnetic field effect that is weak 
enough not to directly affect the superconducting order 
parameter and put aside the field coupling with conduc- 
tion electrons. In the present case, the Zeeman splitting 
of the impurity triplet states reduces the effective pairing 
interaction and the a u -t u scattering depends on the field 
direction. Consequently, the reduced T c exhibits the field 
orientation dependence. 

It is convenient to choose the quantization axis in the 
direction of the applied magnetic field H. The Zeeman 
splitting is expressed as 

{±\-J-h\±) = Th t , (0| — J • h\0) = 0, (20) 

where h t = (5/2)/i and h = \h\ (h = gjp^H: gj is the 
Lande g factor) for the T 5 triplet, and \n) (n = +, 0, — ) is 
an eigenstate oi—Jh expressed by a combination of the 
three T 5 states in eq. (2). For the interchange between the 
singlet and triplet states, the T n operators in eq. (5) cor- 
respondingly depend on the field directions (h x , h y , h z ), 
as indicated in eqs. (A-ll) and (A-12). The magnetic field 
effect only modifies /^= w ,a in the gap equation [eq. (13)] 
as 

ft -+ C h fs(x ) + 1(1 - C h )[fs(x + ) + f 6 (x-)}, (21) 

where the energy difference x n (n = +, 0, — for the ex- 
cited triplet) and the anisotropy coefficient Cy L are de- 
fined as 

Sn - S Fl 



Xri 



2T C 



c h = ^ (rii/sioxoi/siri) 



(22) 
(23) 



respectively. Here, Trf „ indicates the trace of [• • • ] in 
calculating with the 4x4 matrices in eqs. (A-7) and 
(A-8). 

4. Application 

In this section, we consider a case of the single a u 
and threefold degenerate t u bands as an extension of the 
above two-band case (see Fig. 2). The inclusion of the 
threefold degeneracy is required in the minimum model 
for the conservation of the cubic symmetry. We assume 
that these bands are combined with each other only 
through the impurity interband scattering effect. Since 
it is shown that T c is not markedly field-orientation- 
dependent under the cubic symmetry, we introduce a 
uniaxial anisotropy effect phenomenologically to check 
how T c is affected by the orbital (multipole) anisotropy 
of impurity scattering. Finally, we give a comment on 
Hc2 {T) in the vicinity of T c that reflects the anisotropic 
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t t t < >a 

u,x u,y w,z ^ ^ u 



Zeeman 
splitting 



I 



multipolar 
coupling 



Fig. 2. Sketch of the inelastic multipolar exchange scattering in 
the a u and threefold degenerate t u band case. Each local t u ,a 
electron participates in the corresponding (+, a) band (a = 
x,y,z), while the local a u electron is transferred to the single 
(— ) band. 



scattering. 

4-1 a u and threefold degenerate t u band case 

When a magnetic field is applied, it is generally ob- 
servable that under the cubic symmetry, T c decreases as 

T c (h = 0) - T c = a 2 h 2 t + a 4 h*, (24) 

where h t — (5/2)h is the Zeeman splitting, a 2 is a pos- 
itive constant, and 0,4 is the field-orientation-dependent 
parameter 



a 4 oc const. + (h* + + h*). 



(25) 



Here, we simplify the mixing of the local t u (x, y, z) 
orbitals and the + band. Instead, for the conservation 
of the cubic symmetry, we introduce the threefold de- 
generate i u -dominant bands denoted by (+,a), where 
a = x,y, z, which hybridize with the /-orbitals most 
strongly in the directions of three principal axes. The 
corresponding a u -t u scattering anisotropy in eq. (A- 15) 
is given by 





= C h (t 


2' 


Ch,y 






= C h (t 


>=-, 
2' 


Ch,z 


= C h (( 


? = 0) 



= 0) 



2' 



1 

6 " 
1 

" 6 



(26) 



1 1 T2 

6 + 2^ 

When the single a u and threefold degenerate t u bands 
are taken into account, the gap equation [eq. (13)] is 
extended to the following eigenvalue problem: 



T, 



lOf 



c() 



T 

± c 



( A+,x ^ 




f A +,* \ 


A+, a 


= A 


A+, y 
A+, z 


V A - J 







F 



A = a s 















1 U),Z 

F A .z 



Fa.x 
Fa,z 



(27) 



Here, A +iQ , (a = x, y, z) is the order parameter for the 
(+, a) band and no interband pairing is taken into ac- 
count. The magnetic scattering strength as is given by 



ni mp N Q J s 

In the presence of a magnetic field, the matrix elements 
of A for £ (= u), A) are given in the same manner as 
eq. (21): 

Ft, a = Ch^Xo) + i(i - C H , a )\fi{x + ) + fax-)]. 

(29) 

The excited triplet energy levels are split as 

<*o - *ri ( h 

, x ± =x oT x h (xh = W 



x 



(30) 



owing to the Zeeman splitting h t = (5/2)h of the 
triplet. When the field direction is so rotated as to pass 
through the [111] axis, both x and y components are 
taken to be equivalent (h x = h y ). In Appendix C, we 
show that the field-direction (ft, z )-dependent term is ex- 
tracted from Fc >a . 

The highest eigenvalue of A determines T c . When a 
magnetic field is absent (x^ = 0), Ft a = fi,a = ff hi 
eq. (B-l) leads to 

log = a s (2f u + Jf* + 3fi) • (31) 

The xo dependence of this equation is similar to f(x) 
in Fig. B-l that corresponds to the two-band case. The 
calculated T c is plotted for various values of the singlet- 
triplet level splitting (<5r 5 — Stx)/T c q in Fig- 3. The 
as dependence is not sensibly affected by the large 
change in crystal field level for (<5p 5 — <5ri) ^ 10T c o- 
The monotonic increase in T c with as oc ni mp explains 
well the Pr concentration x dependence of T c measured 
in (Lai_ 2; Pr :! ;)Os4Sbi2 at the relatively small x values 
where the Pr ions are regarded as impurities. 15 ^ 

In the presence of a magnetic field, the effect of 77^ in 
eq. (C-2) appears at the field-orientation-dependent T c in 
eq. (24) such that 



0,4 



(32) 



(33) 



where h x = h y is retained. This leads to 

°4[001] — a 4[lll] _ ^ 
a 4[110] — &4[lll] 

In the present model, we find that the difference (TJ^ooi] — 
F c [m]) is of the order of /i 4 , less than 10 -4 T C even for 
a larger magnetic field Xh > 1 since 1 0,4 1 is extremely 
smaller than a 2 . Therefore, r/j is negligible and T c can be 
regarded as isotropic, which is determined by 



T c T c 
^--log — 

J c0 i c0 



a s 2F L 



(34) 



with the same form as eq. (31). The decrease in T c with 
increasing h comes from 

F(x , Xh) - /(xo) = ^f"(xo)xl < 0, (35) 
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Fig. 3. T c /T c q of the s-|--wave state as a function of the inter- 
band magnetic scattering strength as for various values of the 
level splitting (<5r 5 — St-l)/Tco = 5 (solid line), 10 (dashed line), 
15 (one-dashed and one-dotted line), 20 (one-dashed and two- 
dotted line), and 25 (two-dashed and one-dotted line). The same 
dependence is obtained for the interband nonmagnetic scattering 
in the s+-|_-wave state. 



where F = —F& + F u and / = —/a + fu > 0, calculated 
using eq. (C-l). It should be noted that T c could be in- 
creased by applying a magnetic field if f(x) had some 
features to satisfy f"(x ) > 0. 

4-2 Uniaxial anisotropy effect 

As discussed above, no distinct feature is found in the 
field orientation dependence of T c in the rotation of the 
field direction. This is due to the equivalency of x, y, 
and z components in the orbital symmetry and the con- 
stant a 2 in eq. (24). Here, we check how the T c devia- 
tion is increased by lowering the crystal field symmetry 
from Oh- To show this explicitly, we introduce a uniaxi- 
ally anisotropic deviation from Oh- Consequently, a<i has 
a linear dependence on h 2 ,, which leads to the field an- 
gle dependence as cos 20 h (cos8h = h z ), reflecting the 
twofold symmetry. 

Here, we represent the uniaxial anisotropy by the sin- 
gle phenomenological parameter v in the gap equation 
[eq. (27)], modifying F^ iCt as 

F i<a (1 - v)F i<a (a = x,y), F e>z -> (1 + 2v)F^. 

(36) 

For the impurity states, only the Zeeman splitting Xh 
is considered. We assume that the main contribution to 
v comes from the change in local orbital hybridization 
amplitude owing to the lowering of the symmetry. If the 
limits v — > 1 and as — > as/3 are considered, eq. (27) is 
reduced to the two-band case in eq. (13). Around the z- 
axis, we find that 02 in eq. (24) is independent of (h x , h v ) 
since the fourfold symmetry is conserved. In the following 
argument, we calculate 



a 2 ( 



Oh) = 2 ( a 2[ooi] + a2[no]) + 2 ( a 2[ooi] - a 2 [no]) cos 26 h 

(37) 

for various v values, keeping h x — h y in eq. (27). As it 
is expected, the T c deviation becomes more distinct as 




Fig. 4. Anisotropy v dependence of a^ifih) for as = 1 and 
(So — Svf)/T C Q = 10. The solid and dashed lines are the plots 
for [001] (0 h = 0) and fe||[110] (0 h = tt/2), respectively. a 2 
exhibits the cos26^ oscillation between 02(0) and 02(tt/2). 



the crystal field anisotropy v increases. Figure 4 shows 
how d2 increases with \v\ for both ft.||[001] (6h — 0) and 
ft.||[110] {Oh = tt/2). The maximum of T c (minimum of 
02) is given by T" c [ooi] when v > and by T c [no] when v < 
0. It is necessary for the larger amplitude of oscillation 
in eq. (37) to introduce a stronger anisotropy For v = 
0.5, the difference (fflamoi — a 2[ooi]) is estimated as ~ 
10 _3 /T C in Fig. 4. This indicates that the amplitude of 
T c oscillation (oc h 2 ) is about 0.1% of T c at h t /T c ~ 1 
and 1% at h%/T c ~ 3 (ht represents the Zeeman splitting 
of the excited triplet). If the crystal field anisotropy is 
much smaller, it is more difficult to obtain the field angle 
Oh dependence of T c . We note that 77^ in eq. (C-2) is 
responsible for the Oh dependence, so that T c at h z = 
l/v3 for v ^ (T c [ m ] in the cubic symmetry for v = 0) 
directly reflects the crystal field anisotropy (i„-orbital 
anisotropy) since 775 vanishes at h z = l/y/3. Therefore, 
the anisotropy of multipolar scattering itself appears in 
the T c deviation from T c (h z — 1/V3). 

In a real system, this anisotropy effect can be observed 
as the field orientation dependence of the upper critical 
field hc2 near T c at h = 0. We can estimate h C 2 in the 
framework of the conventional Ginzburg-Landau (GL) 
theory, taking into account a diamagnetic effect in the 
GL expansion of a free energy. 32 ) h C 2 is given by At (= 
1 - T/T c ) expansion as 33 ) 



h c2 = Cl At + c 2 (At) 2 



(38) 



with ci > 0. In the present case, the effective pairing in- 
teraction depends on the crystal field energy levels cou- 
pled to the magnetic field, and T c is reduced as 

T c = T c {0)~a2(0h)h 2 C 2 [T C (Q) = T c (h = 0)] . (39) 

Substituting it for T c in eq. (38), we obtain 

T 



Cl 



1 



T c (0) 



C2 



r 3 
c l 



T c (0) 



a 2 (6 h ) 



T 



T c (0) 
(40) 



When the cubic symmetry is conserved, h c2 is invariant 
against the rotation of the field direction since a2(0h) is 
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a constant for v = 0. If a uniaxial anisotropy is applied, 
for instance, by pressure measurement, the field-angle- 
dependent h C 2(T) lines could be observed in the super- 
conducting phase. This is more promising for a large ci 
that is related to the GL parameter. This provides us 
with conclusive evidence of the multiband picture pro- 
posed in our theory The multipole moments in the /- 
electron states play an important role in the anisotropic 
h c2 {T), which is analogous to the quadrupole ordering 
transition temperature Tq(H) with the field orientation 
dependence, as observed in PrPb 3 34 ' 35 ) or as predicted 
for CeBg in a high magnetic field. 36 ) 

5. Conclusion 

We have studied a magnetic field effect on T c increased 
by inelastic scattering for the singlet-triplet configuration 
in impurities. The T c increase can be expected for either 
the interband magnetic scattering in the s±-wave state 
or the nonmagnetic scattering in the s++-wave state. We 
have focused on the r 5 -type octupolar exchange scatter- 
ing in the former. The same argument is applicable to 
the F5-type quadrupolar exchange scattering in the lat- 
ter. Owing to the anisotropy of multipolar scattering [r]£ 
in eq. (C-2)], T c exhibits the field orientation dependence. 
Since the T c deviation (oc h 4 ) is very small in the cubic 
symmetric environment, we have introduced a uniaxial 
anisotropy into the degenerate t u orbitals to clarify the 
multipolar scattering effect on T c . The multipolar scat- 
tering anisotropy itself appears in the twofold symmetric 
oscillation of T c with the rotation of the field direction, 
and the amplitude of oscillation increases proportionally 
to h 2 . This can be confirmed by observing the splitting 
of an H c2 (T) line near T c with a change in field angle, 
although a uniaxial anisotropic field is required to pro- 
duce a clear difference between the two field directions, 
as discussed in §4.2. 

Thus, we have clarified the roles of multipole degrees of 
freedom in multiband superconductors. The key is that 
the local orbital exchange scattering is directly connected 
to the interband scattering via the hybridization of f- 
orbitals with each band. The local orbital anisotropy 
in eq. (C-2) appears clearly at the field-orientation- 
dependent T c , which is related to the field orientation 
dependence in eq. (A-15) for the octupolar scattering. 
This is not expected for the spin exchange scattering in 
a single band. If each t n (rj — z, ±) is replaced by the 
Pauli matrix a v (a± — a x ±i<7 y ) in eq. (A-14), the calcu- 
lated Ch becomes a constant, indicating that the spin ex- 
change scattering is invariant against the rotation of the 
field direction. Therefore, the field-orientation-dependent 
T c , i.e., the anisotropic H c2 , would be positive evidence 
for the multiband picture proposed here. 

In a real system, the quadrupolar scattering coexists 
with the octupolar scattering for the Oh Ti singlet and 
r§ triplet configurations. 7 '' The quadrupolar scattering 
strength oq is defined as as- Unlike the magnetic case, 
the interband nonmagnetic scattering suppresses T c in 
the s±-wave state, which is due to the sign reversal of 
Fa,u (a = x,y,z) in the gap equation [eq. (27)]. The 
competition of magnetic and nonmagnetic scattering ef- 



fects can be taken into account by modifying 

asF U!a -> (a s + a Q )F Ut0l , a s F A , Q -> (a s - a Q )F At0l 

(41) 

in eq. (27). The calculated T c decreases linearly with 
OiQ/ois in the s±-wave state (A +;Ct A_ < 0). For the 
strong spin-orbit coupling, aq/as = 1/9 is derived from 
the Anderson model including the effective exchange 
scattering of the J — 5/2 (r7©rg for Oh in Appendix A) 
electrons due to a single impurity. 7 ' 31 ) This ratio of scat- 
tering strengths leads to only an approximately 10% re- 
duction in T c . Therefore, the T c increase holds for the 
dominant magnetic scattering in the s±-wave state. On 
the other hand, in the s ++ -wave state (A +;Ct A_ > 0), 
the nonmagnetic scattering contributes to the T c increase 
for uq >• as- 7 ^ Experimentally, it has not been estab- 
lished which scattering is more dominant, magnetic or 
nonmagnetic, as the Pr impurity effect on the LaOs4Sbi2 
superconductor. We would like to point out that T c can 
also be suppressed by intraband magnetic scattering that 
corresponds to the spin exchange scattering in the single 
a u band that we have neglected. Since the a u electrons 
are coupled only to the excited triplet in the strong spin- 
orbit coupling case, 31 ) the intraband scattering effect on 
T c is negligibly small. 

As mentioned in §1, the Th symmetry is another fea- 
ture of the skutterudites, leading to the mixing of the Oh 
T 4 and T 5 wave functions of the triplet states as 19, 20 ) 

|ri 2) ) = v / i-rf 2 |r 5 ) + d|r 4 ), (42) 

where d represents the deviation from the Oh symme- 
try. Since d is relatively small, the Th effect modifies ry^ 
slightly in eq. (C-2). On the other hand, the magnetic 
field couples the ground state and one of the triplet states 
via the Van Vleck process, which shifts the energy dif- 
ference x n in eq. (22) by ~ d 2 h 2 /(5o — (5pJ. Thus, the 
Th deviations from Oh only give minor corrections to the 
present results as long as the magnetic field is not large. 
It should be noted that the magnetic field effect on T c is 
more sensitive to the local hybridization of /-electrons 
with conduction bands that we have assumed to be the 
strongest in the directions of three principal axes, leading 
to the maximum or minimum T c in the field orientation 
dependence. 

Finally, we would like to refer to a few experimen- 
tal studies to elucidate the crucial roles of the localized 
Pr 4/-electrons in the superconductivity. In Pr a; Os4Sbi2 
synthesized under a high pressure, the resistivity and 
magnetization data show the close correlation between 
the Pr singlet-triplet energy splitting and Tc- 37 - 1 The re- 
cent nuclear magnetic resonance study indicates the rele- 
vance of magnetic multipole (dipole and octupole) fluctu- 
ations for mass enhancement in PrOs 4 Sbi 2 . 38 * 1 For com- 
parison with our scenario in the future, systematic ex- 
perimental studies of H C 2(T) in (Lai_ a; Pr a ;)Os4Sbi2 are 
highly desired, following the detailed analysis of H c2 (T) 
in PrOs4Sbi2 reported previously. 39 ^ The observation of 
the anisotropic H c2 (T) is worth testing under the uniax- 
ial pressure that could enhance the Pr multipolar scat- 
tering anisotropy in a magnetic field. 
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Appendix A: Interband Impurity Scattering 
with a u -t u Orbital Exchange 

The single /-electron states with the J = 5/2 (J z = 
5/2,3/2, •• • ,—5/2) total angular momentum are classi- 
fied into the Oh symmetric states as 40 ' 



;,3/2 



13/2) 



5/2), 



|r 8 ,i/ 2 > = |l/2), 

|r 8 ,-i/ 2 > = -I - 1/2), 

|r 8 ,_ 3/2 > = - 3/2) + A /||5/2>, 



l r 7,l/2> - 
l r 7 -1/2) = 



(AT) 



-3/2)-W-|5/2) 



>2) 



(A-2) 



gl-5/2). 



The fourfold degenerate r 8 wave functions mix with the 
threefold degenerate t u (x,y,z) orbitals as 

1 



r 8 ,-3/2> ^=(|a;,4.>-%,4.», 



'8,3/2) --^=(|s,t>+%,t», 



l r 8,l/2) ^ 



V3 



- 1/2> " Ts 



V2\zA)-^(\x,l)+i\y,U) 
V2\z,i) + ^=(\xA)-i\yA)) 



(A-3) 



In the same manner, the doubly degenerate electrons 
are directly transferred to the single a u (xyz) orbital as 



|r 7 ,i/2> i \xyz,V, 
|r 7 ,-i/ 2 > i \xyz,i). 



(AT) 



In eq. (5), the T v octupole operators couple both the 
Ti and T 5 states, and their matrix expressions are given 
by 7 ) 



(A5) 



T 1 — 



( 





1 


\ 


























1 

















V 








J 












/ 











1 \ 


rl = 


V2 




-1 










































J 



(A6) 



The interchange of the singlet ground and triplet excited 
states occurs via exchange in and T 8 local electrons 
hybridizing with the a u and t u bands, respectively. In- 



troducing ip = (V'+t ip-t ip-lY f° r the electrons, 
where W<7 is the field operator for the /j = + (t u ) and 
(i = — (a u ) bands with the spin a (=|,|), we present the 
typical octupolar exchange operators with the following 
matrix expressions: 7 -* 



: 2 






c 

\ 8*6* 



.1 
1 — 

2 






— c 


— s 2 e 





s 2 e* 


— c 


-s 2 e 


-i0 





c 














— y/3s\e 











s 2 e* 


c 








v/3sie^ 






(A-7) 



c 

-« 2 e i * 



(A- 



In the matrix elements, c, si, and s 2 are given by 



. = \j - cos( 



si = — p sin#, 
V2 



s 2 = —= sin#, 
v6 



(A9) 



where 9 takes arbitrary values as well as 0, which are 
related to the mixing of the local t u (x,y, z) symmetric 
orbitals and the + band shown as 



(x\+) : (y\+) : (z\+) 



sin V cos < 



sin v sm c 



: cos 9. 
(AT0) 



Here, |+) represents a partial wave of the + band elec- 
trons with the t u symmetry at an impurity site. The 
details of impurity scattering due to such a multipolc 
as the octupole are described for the f 2 singlet-triplet 
configuration in the previous paper. 7 ) 

In the presence of a magnetic field, the T n operators 
are expressed as 



T x = 



T, = 




b 

(a+ - a-) 
-b 



-V2a+e^ 
V^be^ 



(a + — a_ 










) 



b \ 







V2be h 






(ATI) 







(A- 12) 



on the basis of (ri,+,0, — ), where (+,0,—) denote the 
excited triplet states that satisfy eq. (20). The parame- 
ters in each matrix element are magnetic-ficld-dcpendent 
and are defined as 

a± = ^(l±ft z ), b = + h\, tan^> h = ^, 

(AT3) 

where h/h = (h x ,h y ,h z ). Equations (ATI) and (AT2) 
for h z = 1 correspond to eqs. (A-5) and (A-6), respec- 
tively. Using the matrices in eqs. (A-7) and (A-8) to cal- 
culate 



<ri|/ s |0> = (a+ 



a-)i z + -^be^H- + -^be- 
y/2 V2 



(AT4) 
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in eq. (23), we obtain 
Ch = \\ + a + a - 



sin 2 6 



- - 2a+a_ 



cos 2 6> 



+ 2^fa+aZ(a + — a_ ) sin 9 cos 9 cos(c/)h 
+ (a + a_ ) sin 2 9 cos 2((f>h — 4>) i 



(A-15) 



which depends on the details of the hybridization of local 
/-electron states with the + band represented by 9 and 
4> in eq. (A- 10). We note the following equality: 



]T Tr fa [<ri|/,s|n)<n|Js|ri> 

n=0,± 



1. 



(A-16) 



To compare the magnetic and nonmagnetic exchange 
properties, we also present the IVtype quadrupole oper- 
ators 

\ 





o / 
/ o 



Q z = 



( o 


i 









(A-17) 



Q+ = Q 



t 



-V2 



1 



V o 










1 \ 




o J 



(A-18) 



which correspond to the T v octupole operators in 
cqs. (A-5) and (A-6), respectively, and the quadrupolar 
exchange operators 









c 


s 2 e 








-8*6* 


c 


c 


-s 2 e"^ 








s 2 e i * 


c 









i— 

2 













c 









(A-19) 

-% ^ 

s 2 e^ 


J 

(A-20) 



which correspond to the octupolar exchange operators 
(i, q ) in eqs. (A-7) and (A-8), respectively. 7 ) Let us take 
si = s 2 = to simplify i z and q z as 

/ o 



tz 




i 

v o 
/ o 



1 

V o 



-i 





1 



\ 

— i 





J 



o \ 
1 


o J 



(A-21) 



They express the spin-independent scattering with only 
orbital exchange. In the band space, it is clear that the 
former f 2 type is magnetic and the latter fi type is non- 
magnetic. 7 ' It is easy to check whether the effective par- 
ing interaction mediated by such impurity scattering is 
attractive or repulsive in the s±-wave state as follows. 



^0.2 - 




Fig. B-l. Plot of f(x). 



Since the s±-wave state is expressed by Af3/5 2 (J 2 in the 
t® a space, the f 2 -type magnetic scattering satisfies 



f 2 (Af 3 /S 2 <j 2 )f 2 = -Af 3 /5 2 <r 2 , 



(A-22) 



which is used for deriving /a in eq. (13). The sign re- 
versal of A indicates that the impurity-mediated pairing 
interaction is attractive. On the contrary, it becomes re- 
pulsive for the fi-type nonmagnetic scattering since 



fi/3 3 (Af3p 2 <7 2 )fi/5 3 = Af 3/ 3 2 <T 2 



(A-23) 



results in the absence of sign reversal of A owing to p 3 
in the particle-hole space that accompanies the nonmag- 
netic scattering. In a similar analysis, one can confirm 
that the nonmagnetic (magnetic) scattering leads to the 
attractive (repulsive) interaction effectively for the s ++ - 
wave pairing. Therefore, T c can be increased by either 
the magnetic scattering in the s±-wave state or the non- 
magnetic scattering in the s ++ -wave state in the case 
of interband impurity scattering with the a u -t u orbital 
exchange. 

Appendix B: Function as the Impurity Effect on 

Each gap equation presented in this paper gives the 
level-splitting x-dependent T c determined by f(x) = 
—fA(%) + fu(x) that is defined in eq. (16). For the calcu- 
lation of f(x) plotted in Fig. B-l, the necessary formulas 
are arranged here: 7,8 ) 

f A (x) = — +A(x)--B(x), 

9 1 

f u (x) = -l + tanhz- -B(x), (B-l) 

where A(x) = Si (x) tanh x and B{x) = 5 2 (x)tanhx are 
derived from the following equations as 

, oo ipfl + n-i—^-ipl- 
s i( x ) = — Re V 

7T *■ — ' 

n=0 









1 




















" 2 
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^)4lmE , — , 2 W , (B-2) 



n=0 



1 .a; 

Z 7T 



and ^ represents the digamma function. 

Appendix C: Multipolar Scattering Anisotropy 
in a Magnetic Field 

It is convenient to divide i<£ iQ in eq. (29) into the h z - 
indepcndent term 



F 6 (x ,x h ) = - [ft(x ) + fe(x + ) + fc(x-)] 



and the /^-dependent term 
1 (\ 



(C-l) 

n 

h 2 r )[2ft(xo)-fr(x+)-ft{x-)], 



so that F^ yCl is rewritten as 



F 



(C-2) 

(C-3) 
(C-4) 

The z component of the field h z (0 < h z < 1) represents 
afield direction: [110] (/i 2 = 0), [001] (h* = 1), and [111] 
(h z = l/\/3). One can see that T c has a field orientation 
h z dependence, which comes from eq. (C-2). For Xh <C 1, 
it is reduced to 



Ft, z = F^(xo,Xh) + 2i] i (xo,Xh)- 



T]£(xo,Xh 



h 2 z tt'(x )xl 



where /" indicates the second derivative. 



(C-5) 
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